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We conjecture (analytically) and demonstrate (numerically) the existence of a fine-structure above 
the power-law behavior of the mass of black-holes that form in gravitational collapse of spherical 
massless scalar field jjj . The fine-structure is a periodic function of the critical-separation (p — p*). 
We predict that the period vj is universal and that it depends on the previous universal parameters, 
the critical exponent, ft, and the echoing period A as zu = A//3. 
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The gravitational collapse of a spherically-symmetric massless scalar-field has two possible end states. Either the 
scalar field dissipates away leaving a flat spacetime or a black-hole forms. Numerical simulations of this model problem 
, Q have revealed an unexpected critical behavior when the initial conditions are close to a critical case p — p* (p 
^ \ is some parameter which characterizes the strength of the initial scalar field, and p* is the threshold value). More 
■ precisely, Choptuik has found a power-law dependence of the black- hole mass on critical separation (p — p*) of the 
"sj" ' form: oc (p — p*) 13 for p > p* , and a discrete echoing with a period, A, (a discrete self-similar behavior) for 

p = p*. 

Subsequently similar critical behavior have been observed for other collapsing fields: axisymmetric gravitational 
wave packets || , spherically symmetric radiative fluids Q| and complex scalar fields [^) . In all these model problems 
the critical exponent (3 turned out to be close to the value originally found by Choptuik 0.37, suggesting a universal 
! behavior. However, Maison [|| has shown that for fluid collapse models with an equation of state given by p — kp the 
\^ • critical exponent strongly depends on the parameter k. 

In this work we conjecture the existance of a small periodic correction, \&[(ln(p— p*)], to the power-law dependence of 
^sO . the black-hole mass. *S> is periodic and its period, w, is universal and it depends on the previous universal parameters 
1 as w = A/ p. Our analytical argument predicts the existence of the fine-structure periodic term and its expected 
, period. The argument is based upon the final stage of a super critical evolution: from the moment when the deviation 
'■ from the exact self-similar critical evolution becomes larger than some given value (and the evolution is no longer 
| | self- similar) up to the horizon formation. We then provide a numerical evidence that verifies the existence of the 
W)' conjectured periodic term, the universality of its period and the relation w = A//3. 

Our numerical formalism is based on the characteristic scheme of Goldwirth and Piran |7j to which we have added 
• . an expansion near the origin which is essential to achieve the extremely high accuracy needed for these computations. 
p\ ' The evolution equations, our algorithm and numerical methods, and the discretization and error analysis are all 
^ , described in a previous paper ^ , and will not be repeated here. 



II. THEORETICAL PREDICTIONS VS. NUMERICAL RESULTS. 



We consider the spherical collapse of a massless scalar field. Choptuik has shown that for a critical parameter 
p* there is a critical solution which has an infinite discrete self-similar behavior. The critical solution by itself does 
not yield the black-hole mass scaling relation in which we are interested. We perturb, therefore, the critical initial 
conditions. This leads to a dynamical instability - a growing deviation from the critical evolution toward either 
sub-critical dissipation or supercritical black-hole formation. 

Let f(u) be a function of u, the time coordinate of an observer at rest at the origin, that characterize the solution 
along the outgoing null geodesic that leaves the origin at u. The function, /, could be, for example, the maximal 
value of M(r, u)/r along this geodesic. Following Evans and Coleman Q and Maison [|| we describe the run-away of 
the perturbed solution from the critical evolution, (described by f c ) as a power-law: 

f(u) - f c (u) = X(u* - u)- a , (1) 

where the critical solution reaches the zero-mass singularity at u — u* . The prefactor A satisfies: A oc (p — p*). 

We assume that the range of validity of the perturbation theory is restricted to some maximal deviation, x, from the 
exact critical evolution, i.e. the evolution is approximately self-similar until / — f c = x- From here on the evolution 
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is outside the scope of the perturbation theory - there is sub-critical dissipation of the field or supercritical black-hole 
formation. In cither case, the evolution from this stage onwards looses its self-similar character. We chose now p > p* 
so that the perturbed initial conditions develop into a black- hole. The time u x {p — p*) required in order to reach the 
maximal deviation is given simply by the relation: 

X(u* u x y a = X . (2) 

Of course, a larger initial perturbation requires a shorter time to reach this value. We define now the logarithmic 
time, T = — ln[(u* — u)/u*], in which the critical solution is periodic. The logarithmic time, T x , which corresponds 
to the loss of self-similarity, is given by 

T x = -a- 1 ln(p-p*) + b k , (3) 

where depends on Xi u * an d k. The index k denotes the family of initial conditions considered. We conjecture now 
that the logarithmic time until the horizon formation, equals to T x plus a periodic term F[(\n(p — p*)]: 

T bh = -a- 1 ki(p-p*) + b k +F[ln(p-p*)} . (4) 
The period, w, is universal and it depends on the previous universal parameters according to: 

m = aA . (5) 

Consider now two different initial-conditions, which lead to n and n+1 echoes respectively (until the deviation from 
the critical self-similar evolution reaches x). These solutions are related to each other by an exact scale transformation 
with a factor e A . The final stages of these two supercritical evolutions, from the stage when the deviation from the 
exact critical evolution reaches x, (and the evolution ceases to be periodic in T), up to the horizon formation, are 
equal up to a scaling transformation. The periodic nature of the function F arises from this final stage: The period 
of the function F is the amount that should be added to the quantity ln(p — p*), in order to reduce the number of 
echoes by one. This will reduce Tbh by A. From Eq. |^ this amount to a period vo — Aa in F. 

This conjecture is verified by numerical simulations of four families of initial data (two neutral and two charged). 
In all those families we have found that —Tht as a function of ln(p — p*) was well fit by a straight line with a slope 
1/a « 0.37. On top of this straight line there was a small modulation. The deviation from a straight line is shown 
in Fig. 1. which provides a numerical evidence for the existence of the periodic term F in Eq. ^. We see that the 
function F is indeed periodic, with a universal period <p = Aa w 4.6. 
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FIG. 1. Illustration of the conjectured universal periodic fine- structure of —Tbh- The quantity [—Tbh — {—Tbh}] is plotted as 
a function of ln(a) , where a = (p — p*) /p* , for the four families. The curves were shifted horizontally (but not vertically) 
in order to overlap the first oscilation of each family with the first one of family (a). (Tbh) is the value of Tbh determined 
from a straight line approximation, i.e. {Tbh) = Const + f3\n(p — p*). The numerical results agree with the predicted relation 
■co = aA » 4.6. 
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We have to relate, now the exponent a to the critical exponent (3 (which describes the power-law dependence of 
the black-hole mass) and then generalize Choptuik's scaling-relation by proving that one should also add a periodic 
term to Choptuik's mass scaling relation. We write T b h in the form: 

T bh = T init +nA + F , (6) 

where Ti n n is the initial logarithmic time required for the system to settles down to a periodic behavior in T, and n 
is the number of echoes. We assume that Tina is independent of (p — p*). Using Eq. |] we obtain: 

nA = -a- 1 ln(p-p*)+4 , (7) 

where dk is a family-dependent constant. We define AfW as the mass after n echoes (note that this is not the final 
black hole mass). Since M decreases in each echo by a factor e~ A we have, using Eq. [?]: 

M< n > = M^e- nA = M^e~ d " {p - p*f , (8) 

from which it follows that (3 = 1/a. 

To obtain, Mbh, the final black-hole mass one should multiply by a periodic function G[ln(p — p*)] which 

measures the change of mass, from the stage when the evolution is no longer periodic in T, until the horizon forms. 
The function G depends only on the field configuration at the moment when the deviation from the exact self-similar 
evolution reaches x (and the evolution is no longer self-similar). Thus, G is expected to have the same value each 
time the system completes another echo, i.e. each time n increases by unity. Using Eq. (6) we find that the function 
G[ln(p — p*)] is expected to have a period of w — A/ (3. Thus, we obtain 

ln{M bh ) = 0hx{p-p*) +c k + ¥pn(p-p*)] , (9) 

where Ck is a family-dependent constant and ^[ln^ — p*)] is a periodic function with a universal period, w. 

Fig. 2 depicts this periodic fine structure for our four families of solutions mentioned earlier. In all four families we 
obtain the basic power law behavior with /3 « 0.37. Fig. 2 displays the deviation of ln(M m h) from this straight line 
as a function of ln(p — p*). The agreement between the four families shows that the fine structure is indeed universal 
with the expected period. 
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FIG. 2. Illustration of the conjectured universal periodic fine- structure generalization of Choptuik's mass-scaling relation. 
ln(m) — (ln(m)) is plotted as a function of ln(a) for the four families, where m = Mt,h/Mi,h,c is the normalized black hole 
mass in units of the initial mass in the critical solution Mf,h,c- (ln(m)) is the value o/ln(m) determined from a straight line 
approximation. The curves were shifted horizontally (but not vertically) in order to overlap the first oscilation of each familiy 
with the first one of family (a). The numerical results agree with the it predected relation zzi — A/f3 « 4.6. 
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One may worry, of course, whether this fine structure is real or could it arise from some numerical errors. In a 
previous paper |3| we have established the stability and convergence of our code with numerous tests. Still becuase 
of the importance of this issue we demonstrate here the physical character of this fine structure. Fig. 3 depicts the 
deviations of ln(M m &) from a straight line as a function of ln(p — p*) for five different calculations with 100, 200, 400, 
800 and 1600 grid points for the same initial data. The same features were found on all grids even though the grid 
sizes differ by a factor of 16. The five curves overlap and all show the same periodic behaviour. Numerical convergence 
(the 800 and 1600 curves are nearer than the 100 and 200 curves, for example) is clearly seen. 
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FIG. 3. ln(m) — (ln(m)} is plotted as a function of \n(a) for family (a) and for it five different resolution grids with 100,200,400, 
800 and 1600 gridpoints. The five curves overlap and all show the same periodic behavior. 



III. SUMMARY AND CONCLUSIONS. 



We have studied the spherical gravitational collapse of a massless scalar-field, both for the uncharged case and 
for the charged configurations. Our main interest was the supercritical (p > p*) feature of Choptuik's solution, i.e. 
the power-law dependence of the black-hole mass on the critical separation. We have shown the existence of a fine- 
structure above this power-law dependence in the form of a periodic term with a universal period, w. We are not aware 
of such fine-structure periodic term in any other phase-transitions in statistical mechanics. Our periodic term with 
its period strongly depends on the discrete echoing character of the critical solution. This discrete self-similarity has 
been seen only in the works of Choptuik Q , Abrahams and Evans pj (collapse of axisymmetric vacuum gravitational 
field) and in our work concerning the gravitational collapse of a charged (complex) scalar field j^]. Abrahams and 
Evans ||] have found an echoing period of A w 0.3. Thus, using our analytical argument, we conjecture that a careful 
analysis will reveal a periodic fine-structure (to the power-law behavior), with a period of w « 0.8, in the model 
problem of the collapse of axisymmetric gravitational wave packets. 
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